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ABSTRACT 

While the importance of dusty asymptotic giant branch (AGB) stars to galactic chem- 
ical enrichment is widely recognised, a sophisticated understanding of the dust forma- 
tion and wind-driving mechanisms has proven elusive due in part to the difficulty in 
spatially-resolving the dust formation regions themselves. We have observed twenty 
dust-enshrouded AGB stars as part of the Keck Aperture Masking Experiment, re- 
solving all of them in multiple near- infrared bands between 1.5^m and 3.1/im. We 
find 45% of the targets to show measurable elongations that, when correcting for the 
greater distances of the targets, would correspond to significantly asymmetric dust 
shells on par with the well-known cases of IRC +10216 or CIT 6. Using radiative 
transfer models, we find the sublimation temperature of T su b (silicates) = 1130 ± 90 K 
and T su b (amorphous carbon) = 1170±60K, both somewhat lower than expected from 
laboratory measurements and vastly below temperatures inferred from the inner edge 
of YSO disks. The fact that O-rich and C-rich dust types showed the same sublimation 
temperature was surprising as well. For the most optically-thick shells (r2.2^m > 2), 
the temperature profile of the inner dust shell is observed to change substantially, an 
effect we suggest could arise when individual dust clumps become optically-thick at 
the highest mass-loss rates. 

Key words: radiative transfer — instrumentation: interferometers — circumstellar 
matter — stars: AGB stage — stars: dust shells 



1 INTRODUCTION 

One of the most dramatic phases in the life of an inter- 
mediate mass star is the Asymptotic Giant Branch, a rel- 
atively short period where a star loses most of its initial 
mass through a dusty wind. Researchers still do not un- 
derstand all the ingredients necessary for producing the 
high mass-loss rates observed during this stage. The mas- 
sive envelopes ejected during this phase are thought to be 
later illuminated during the planetary nebula stage, a stage 
where most stars show stro ng bipolar circumstellar struc- 
tures (jBalick fc FranklbOOj ). 

Following the advent of infrared detectors, early workers 
made simple spherically-symmetric models of dusty shells 
around large samples of AGB st ars fitting only to the spec- 
tral energy distributions^. g., iRowan-Robinson fc Harris] 
1982, 1983b. a). M-typc stars are typically surrounded by 
dust shells composed of amorphous silicates while C-stars 
have carbonaceous dust. These early workers were able to 



show that dust condensed around 1000 K within a few stel- 
lar radii of the stars and also estimated mass-loss rates typi- 
cally 10~ 6 Mq /yr and as high as 10~ 4 Mq /yr. More recently, 
llvezic fc Elitzurl (|l995l ) developed the code DUSTY to study 
dust shells in a systematic way and made models for a large 
sample of stars, again fitting just the spectral energy distri- 
bution. 

The simple picture of spherically-symmetric and uni- 
form mass-loss was challenged by the observations of the 
Infrared Spatial Interfe rometer (ISI) . a lon g-baseline mid- 
infrared interferometer (|Danchi et aLlll994h . These workers 
found a diversity of shell morphologies with some red gi- 
ants showing episodic dust shells ejections and others with 
a more continuous distribution of dust. A more dynamic 
and asymmetric vision of mass-loss fit into debates into the 
origins of bipolar symmetry in planetary nebulae. High an- 
gular resolution near-infrared speckle and aperture mask- 
ing on 8-m class telescopes were able to image fine de- 
tails on some dust shells, such as the prototype carbon 
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star IRC +10216 (e.g. ITuthill et all l2000al) . An e laborate 
model was presented by Men'shchikov et al.l (|2002f ) arguing 
for complex, spatially-varying dust properties and density 
structures. While IRC +10216 shows complexity within the 
inner few stellar radii, it is unclear if these structures rep- 
resent global asymmetries or just weather conditions of the 
dust formation process observed in situ. 

Here we present the full dataset of dust-enshrouded 
giants observed with the 10-ye ar project called the Keck 
Aperture Masking Experiment (|Tuthill et all 1200001 ). This 
experiment delivered well-calibrated spatial information on 
the scale of ~50 milliarcseconds (mas) in the astronomical 
K band (Ao = 2.2/im), enough to resolve all the dusty tar- 
gets presented here and to measure their dust shell sizes and 
asymmetries. This paper includes 20 objects with observa- 
tions in typically 3 wavelengths ranges, 1.65/im, 2.2/im, and 
3.1/xm. We have also extracted photometry to construct co- 
eval near-IR spectral energy distributions - an important 
factor since these objects pulsate and show large variations 
in flux on yearly timescales. Lastly, we used a radiative 
transfer code to fit each epoch of each target star using 
simultaneously the NIR photometry and multi-wavelength 
angular size information from Keck masking. 

The primary goals of these observations and modeling 
efforts are to measure the physical characteristics of a large 
sample of the most extreme dusty AGB stars, to address the 
question of the onset of circumstellar asymmetries, to deter- 
mine any differences between silicate and carbon-rich dust 
shells, and to constrain the optical properties of the dust 
particles themselves. Lastly, this publication marks the final 
large data release of AGB star data from our diffraction- 
limited Keck masking experiment and we anticipate this 
work will provide a rich dataset for more detailed modelling 
efforts by other workers. 



2 OBSERVATIONS 

2.1 Overview of Observations 

Our observations consist of photometric and visibility data 
taken on 20 different stars at the W.M. Keck observatory 
between December 1997 and July 2002. The wavelengths 
at which these stars were observed and the properties of 
the corresponding filters are listed in Table 1. A listing of 
the observed stars, segregated into carbon-rich and oxygen- 
rich groups, along with their basic properties can be found 
in Table 2. Most stars were measured at more than one 
epoch during this time span allowing for robust internal data 
quality checks. 



2.2 Photometric Data 

Aperture masking procedures consist of alternating target 
and calibrator observations that allow for basic photome- 
try in m ost observing conditions. As part o f the standard 
pipeline (|Monnierjll999l ; ITuthill et al.ll2000bl ) we performed 
aperture photometry on each object, allowing the difference 
in magnitude (Amag) between the target star and calibra- 
tor star to be measured. The Vizier catalog service, most 
often referencing t he Catalogue of Infrared Observations 
jGezari et al.lll999l ) and 2MASS (|Cutrill2003l ). was used to 



determine magnitudes at infrared wavelengths for the cali- 
brators. Interpolation was used between wavelengths found 
in the catalogues and the wavelengths at which our data was 
taken. Occasionally no mid-IR measurements were available 
for some calibrators and we used the calibrator spectral type 
and the K band flux to estimate the flux density at these 
longer wavelengths. 

As a data quality check we compared our photometry 
with 2MASS and found good general agreement, although 
strict agreement was not expected since our targets are 
highly variable and there is some difference in beam sizes. 
We estimated the error on the photometry points at 10% 
based on night-to-night variations. However, there were in- 
stances when we assigned larger errors (between 10 and 32%) 
due to saturation of the 2MASS photometry used for the 
calibrator, intrinsic variability of the calibrator, or effects 
of cirrus clouds in some of the original data. Indeed, there 
were some nights too contaminated by variable clouds to 
allow photometry to be extracted at all. 

Table 3 is a journal of observations, including the ob- 
serving date(s), the filter(s) used, the aperture mask(s) used, 
and calibrator star name. We have compiled the adopted 
calibrator properties in Table 4. 



2.3 Visibility Data 

2. 3. 1 Methodology 

Our group carried out aperture masking interferometry at 
the Keck-1 telescope from 1996 - 2005. We have published 
images and size measurements with (at the time) unprece- 
dented angular resolution on topics ranging from young 
stellar objects, carbon stars, red supergian ts, and photo- 



spheric diameters of Mi ra variables (e.g., Monnier et al.l 



1 19991 ; ITuthill et ai]|2000al lbl; iDanchi et alj|200ll ). 

The NIRC cam era with the image magnifier 
ijMatthews et alj Il996l) was used in conjunction with 
the aperture masking hardware to create fringes at the 
image plane. The data frames were taken in speckle mode 
(Tint=0.14s) to freeze the atmosphere. In the work pre- 
sented here, multiple aperture masks and bandpass filters 
were employed. After fiat-fielding, bad pixel correction, and 
sky-subtraction, Fourier methods were used to extract fringe 
visibilities and closure phases from each frame and averaged 
in groups of 100 frames. Absolute calibration to account 
for the optical transfer function and decoherence from 
atmospheric seeing was performed by interleaving science 
observations with measurements of unresolved calibrators 
stars. At the end of the pipeline, the data products are 
purely interferometric as if obtained with a long-baseline 
interfero meter. A full descripti on of this exper i ment can be 
found in ITuthill et all (|2000bl ) and iMonnierl d 19991). with 
furthe r dis cussion of sys t emati c errors in iMonnier et al.l 
l|2004l ) and lMonnier et al.l (|2007l ). All V 2 and closure phase 
data are available from the authors; all data products 
are stored in the FITS-based, optical interfer ometry data 
excha nge format (OI-FITS), as described in IPauls et al.l 
l|2005f ). 
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2.3.2 Basic Results 

Before undertaking radiative transfer modeling, we provide 
the results of basic geometrical analysis of the visibility 
data. The simplest representation of the data is generally 
a circularly-symmetric Gaussian envelope, a useful model 
to give a characteristic size to the emission. Table 5 pro- 
vides the visibility intercept (Vb) and the Full-width at Half- 
maximum (FWHM) for the best fit for all datasets, including 
the reduced \ 2 - Errors are generally dominated by systemat- 
ics related to the calibration procedure (i.e., seeing variation 
between source and ca librator visits) and w e have used the 
relations established in M onnier et alj (120071 ) to quantify our 
errors. In some cases, there was evidence of two components 
to the visibility curve and we have also fitted a slightly more 
complex model of a point source plus a Gaussian envelope 
to all epochs. Table 6 contains the best fit parameters of 
the 2-component model, including the estimated fraction of 
light in the point source (/ po i n t) and the fraction of light in 
the Gaussian envelope (/Gauss). 

In addition, we fitted each object with a 2-dimensional 
Gaussian function in order to search for signs of asymme- 
try. Objects with observed asymmetry are marked with an 
asterisk in Table 5 . Table 7 lists all the object with con- 
firmed asymmetries and we include the amount of elongation 
( FWHM m ° J °' ) and the position angle (degrees East of North) 
of the major axis. Here we have used the spread of measured 
position angles between wavelength channels and epochs to 
estimate the PA error. We will discuss further these findings 
in H 



3 DUST SHELL MODELING 
3.1 Introduction 

The objects in our study all have spectral energy distribu- 
tions that peak in the infrared. Indeed, these stars are sur- 
rounded by dust shells that absorb the stellar light and then 
reemit the energy in the infrared. In order to extract physical 
characteristics of these dust shells (i.e., optical depths, tem- 
peratures, etc), we must be able to compute how the dust 
will absorb, scatter, and reemit the energy from the star. We 
accomplish this w ith the radiative-transfer model DUSTY 
i|lvezic et al.ll 19991 ). While DUSTY is limited to calculations 
in spherical symmetry, we established in the previous section 
that most of our objects show only mild signs of global asym- 
metries; however, we caution that our results will be suspect 
for the most asymmetric of the targets listed in Table 7. 
Given a small number of input parameters, DUSTY can 
quickly compute synthetic photometry and intensity pro- 
files for dust shells. These outputs can then be compared to 
the data that we have experimentally obtained. 



3.2 Model Description 

We applied a uniform procedure for fitting all of our objects. 
Here we discuss which properties were held fixed and how 
we explored a grid of the key dust shell parameters. 

We begin with the central star. At the beginning of 
our study we used a featureless Planck blackbody spec- 
trum, however we came to realise that a blackbody spec- 
trum is a rather poor approximation for the extremely late- 



type giants in our sample due to strong molecular absorp- 
tion bands. Most notably, the HCN absorption feature of 
carbon-rich stars sits directly at the PAHcs (3.0825 /im) 
wavelength, where we have many observations. Because of 
the severe optical absorption of the dust, spectral types are 
not known for most stars in our sample and we have adopted 
an effective temperature of 2600K for all stars, which is 
as cool as we could find converged synthetic spectra. For 
the carbon stars we used a MARCS model as described in 
iLoidl et all (|200ll ) and for the M-gi ants we used a PHOENIX 
NEXTGEN model as described in lHauschildt et al.l (|l999h . 
The medium-resolution synthetic spectra from these sources 
were smoothed before input into DUSTY. Unfortunately, we 
do not have useful distance estimate to our sources - so we 
adopted a distance of 1000 pc and interstellar reddening 
of Eb-v = 0.5 for all objects. We note that the dust shells 
around the stars absorb nearly all of the energy from the cen- 
tral source, acting as a kind of calorimeter. Thus, while our 
2600K estimate for the central star temperature is crude, we 
expect the bolometric luminosity (for assumed d = 1000 pc) 
to be more accurate. However in practice our luminosity 
estimates are poor due to uncertainties in the dust shell op- 
tical depth and the fact we are not integrating the whole 
observed SED throughout the mid- and far-infrared. 

Based on the shape of the SED (and the presence 
of a silicate feature in IRAS-LRS spectra), we deter- 
mined each star to have either carbon-rich dust or silicate- 
rich dust. Ba sed on this assignment, we chose amor- 
phous carbon ([Hannerl Il988l ) or warm amorphous silicates 
dOssenkopf et al] 1 19921) re spectively in the DUSTY model 
setup. Speck et al. ( 20081 ) discussed how silicates close to 
AGB stars could quickly anneal to crystalline grains but a 
full exploration of optical constants for different grain types 
was beyond the scope of this work. For the grain size distri- 
bution, we adopted the standard MR N power-law grain s ize 
distribution between 0.005-0.25 /im (|Mathis et alj|l977l ); a 
later exploration of larger grain size s did not systemat ically 
improve fits (also see discussion bv lSpeck et al1l2009D . An- 
other property of the dust shell we fixed is that the dust 
density follows a r~ 2 power-law, corresponding to constant 
mass-loss rate. 

Lastly, we come to the parameters of the model that 
are not fixed: the temperature of the dust shell at the inner 
boundary, Td us t, the radius of the star, Rstar, and the K- 
band optical depth T2.2 M m of the dust shell (as integrated 
along the line-of-sight from the observer to the star). In the 
next section, we explain our fitting procedure. 



3.3 Fitting Methodology 

We explored inner dust temperatures Td us t between 400K- 
1500K. This range explored both the high temperatures 
thought to be prohibitive of dust creation and low tem- 
peratures too cool for steady-state dust production. Note 
that when setting up a model in DUSTY, one does not 
specify the inner radius of the dust shell: this quantity is 
calculated based on the luminosity of the star and the spe- 
cific inner shell dust temperature Tdust. In terms of optical 
depth, we explored T2.2 M m between and 9. This range pro- 
vided a full fitting region for our objects and values of T2.2 M m 
much above 9 were too computationally expensive. Finally, 
'Rstar was recognised to simply be a scaling factor for the 
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model outputs and could easily be optimised for every pair 
of (Tdust, T2.2 M m)- Because the DUSTY calculation was fast 
and we only had to optimise over a few parameters, we chose 
to carry out an exhaustive grid calculation over all (Tdust, 

T2.2fj.rn)- 

For each location in the grid we calculated the model 
SED as well as the radial intensity profiles. We calculated 
a x 2 based on both our coeval near-infrared photometry as 
well as Keck masking visibility curves. For the SED, we also 
used including V-band magnitudes in our fit with a very low 
weight to ensure that the optical depths were not too low 
(important especially when for objects without photometry 
in all three near-IR wavelength bands). When calculating 
the % 2 f° r the visibility curves, we adopted the following 
procedure. Because the y-intercept of our observed visibil- 
ity data can fluctuate ±5% due to seeing calibration, we 
normalised each visibility to 1.0 at zero baseline before fit- 
ting. Also, we weighted the visibility points so that the SED 
and the visibility data were separately given equal weight 
in the final reduced x 2 ■ We purposefully chose not to in- 
clude longer wavelength SED measurements, such as IRAS 
data, in our fitting. By fitting only to near-IR photometry 
and near-IR spatial data we can isolate and only probe dust 
emitted within the last few decades. This allows us to keep 
the model as simple as possible and enhances the validity of 
our assumption of constant mass loss rate (ie., p oc r~ 2 ). 

Once the grid calculation over inner dust temperature 
Tdust and T2.2 M m was completed, the x 2 surface was used to 
estimate the best-fit parameters. The uncertainty estimates 
were produced by considering the region where the reduced 
X 2 was less than 2, a highly conservative criterion that re- 
flects the highly-correlated errors in our datasets. In cases 
where the best-fit x 2 is above 1, we scaled the x 2 results by 
the best-fitting value before estimating the parameter uncer- 
tainties. The best-fitting parameters and their uncertainties 
are compiled in Table 8. 

In addition to providing the fitting results in tabulated 
form, we also include here a series of figures which graphi- 
cally represent the new data, modeling results, and the x 2 
surface in our grid. These plots can be found in each of 
Figures 1-20. The first panel in each figure contains the 
observed near-IR photometry and best-fit model SED. The 
second panel in each figure contains the multi-wavelength 
visibility curves averaged azimuthally along with the model 
curves. Finally, the third panel shows the x 2 surface in the 
(Tdust, T2.2fim) plane. We have grouped all the epochs for 
the same object together so one can see the self-consistency 
in the derived dust shell parameters - indeed, consistent 
dust shell properties were recovered when fitting to different 
epochs, despite large changes in the central star luminosity 
due to pulsations. 

One of the most important results to take away from 
these panels that we clearly break the standard degener- 
acy between dust temperature and optical depth. This is 
because of our new spatial information - by measuring the 
sizes of the dust shell at various wavelengths we can simul- 
taneously constrain the temperature and optical depth. In 
the past, one typically had to choose an inner dust tem- 
perature based on physical arguments concerning the dust 
condensation temperatures of various dust species. Here, we 
see that the inner dust temperature can be constrained in- 



dependently from other parameters and the implications are 
discussed further in the next section. 

While the simultaneous fits to the near-IR SED and 
visibility data were generally acceptable, we found the fits 
to the shortest wavelength visibility data at H band were 
systematically worse. Since this band is most sensitive to 
scattering by dust, we explored modified dust distributions, 
especially using larger grains; we did not find systematic 
improvements to the fits by altering dust size distribution 
from MRN or by using other dust constants. 



4 DISCUSSION 

Our survey provides the first constraints on the asymmetry 
of the dust shells for such a large sample of dust-enshrouded 
AGB stars. We found that 4 out of 7 M-stars and 5 of 13 
C-stars showed evidence of dust shell asymmetries, with 
dust shell elongations between 10% and 40%. While this 
level of asymmetry may sound mild, it actually (quanti- 
tatively) compares to the level of asymmetry that would 
be expected for the most asymmetric dust shells known if 
placed at 1 kpc. Fo r instance, we know that IRC + 10216 
i|Tuthill et alj|2000a| ) and CIT 6 (|Monnier et al.ll2000t ) have 
dramatic global asymmetries in their dust shell, detailed 
imaging made possible by virtue of their proximity. If we 
placed these targets farther away, we would not be able to 
image the detail but they would appear ~20% elongated, 
similar to the degree observed here in 45% of our sample. For 
CIT 3 , we confirm t he asymmetries seen b y lHofmann et al.l 
l|200ll) and note that I Vinkovic et all l|2004l showed that the 
20% elongation could be explained by a bipolar outflow . 
That said, clumpy dust formation (|Fleischer et al.1 11992T ) 
might also cause stochastic variations in the inner dust shell 
geometry that could appear as short-lived elongations. Mid- 
infrared observations with long-baseline interferometers (e.g, 
ISI, VLTI-MIDI) should focus on these targets to determine 
the nature of the asymmetries. In addition, long-term moni- 
toring of these dust shells will help settle debates concerning 
when the environments of evolved stars develop large scale 
asymmetries commonly revealed in the later planetary neb- 
ula stage. For instance, a long-term asymmetry in a constant 
position angle (as judged by linear polarization or spatially 
resolved data) would be a sign of a global bipolar mass-loss 
asymmetry and not just weather. 

In order to look at dust shell properties for our full 
sample, we have plotted the inner edge dust temperature 
Tdust vs total dust shell optical depth r2.2pm for all our 
targets. Figure 21 shows these results split into O-rich 
and C-rich dust type. For K-band optical depths below 2, 
we find the sublimation temperature of T su b (silicates) = 
1130±90K and Tsub (amorphous carbon) = 1170±60K, both 
somew hat lower than expected from laboratory measure- 
ments (|Lodders fc Feglevl I1999T I and vastly below temper- 
atures inferred from the i nner edge of YSO d isks (~1800K, 
iTannirkulam et al. | |200S| ; iBenistv et alj|20ich . One compo- 
nent to the observed lower dust temperature could be due 
to the fact that the central star varies in luminosity by 
about a factor of 2 during the pulsation cycle and we see 
the dust cooler than the condensation temperature during 
phases away from maximum light. 

The Tdust vs optical depth T2.2 M m diagram (Figure 21) 
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also shows no statistically-significant difference between O- 
rich and C-rich dust types, counte r to expectation of highe r 
temperatures for carbon-rich dust (|Lodders fc Feglevll 19991 ). 
We recognise that our simple dust shell modeling may not 
lead to accurate estimates of the dust sublimation temper- 
ature if the inner dust formation environment radically de- 
parts from a power law density distribution, perhaps due 
to pulsations, timescale for dust formation, or multiple dust 
species. Interestingly though these concerns would likely af- 
fect C-rich and O-rich shells similarly and so the lack of a 
clear difference in sublimation temperatures between these 
dust types appears robust. 

The other important feature of Fig ure 21, Xdust vs op- 
tical depth T2.2 f im, is the apparent temperature at the inner 
edge of the dust shell gets lower and lower with increasing 
optical depths above 2. This appears true for both C-rich 
and O-rich shells. Here we do not believe we are seeing an 
actual reduction in the dust sublimation temperature, but 
rather a change in the temperature profile in the inner dust 
formation zone due to a breakdown in the assumption of a 
spherically-symmetric r~ 2 density power law. We have am- 
ple evidence that dust formation is cl umpy, as has been im - 
aged in great detail for IRC +10216 (|Tuthill et al.ll2000al ). 
but these clumps have been shown to have a relatively weak 
affect on the temperature structure for low optical depths. 
Next we further explore how a clumpy dusty environment 
could change the temperature profile of the dust shell when 
the individual clumps become themselves optically thick to 
the stellar and even hot dust radiation field. 

Clumpy structures are seen to evolve in 2-D models 
of d ust shells due to se lf-amplifying density perturbations 
(e.g. IWoitke et ai1l2000l ). First optically thick dust regions 
form and these regions cast shadows on the dust behind 
them. Consequently, the temperatures decrease by 100's of 
degrees K and this allows for a higher rate of dust forma- 
tion in these shadow regions. Scattering and re-emission of 
light by the optical thick regions increases the intensity of 
radiation between them and eventually the light escapes 
through the optically thin regions in between the optically 
thick regions. Thereupon, the temperature within the opti- 
cally thin regions increases, which decreases the rate of dust 
production. These processes thus amplify the initial homo- 
geneities until large-scal e clumpy structures start to form, 
such as "dust fingers" l|Woitke fc Niccolinil l2005h . Indeed, 
IWoitke fc Niccolinil (|2005F did see average dust temperatures 
to be reduced due to these opacity effects but at much weaker 
level than we see in Figure 21. Realizing that our data reveal 
a strong effect only a t r's s everal times larger than probed by 
IWoitke fc Niccolinil (120051 ). we suggest that dust shadowing 
effects get dramatically stronger when individual clumps be- 
come optically thick to both stellar radiation as well as hot 
dust emission. A 3D radiative transfer calculation of a dusty 
dust shell could validate or disprove this explanation. 

In conclusion, our large sample of spatially resolved 
dust-enshrouded stars have led to new insights into the late 
stages of AGB star evolution. We find levels of dust shell 
elongations that point to significant asymmetries in nearly 
half of our targets. Our spatial and SED data combined 
has eliminated some model degeneracies, and we now have 
the best constraints on the actual sublimation temperatures 
for dust forming in this outflows, finding lower temperatures 
than expected from terrestrial experiments and not confirm- 



ing the large difference expected between carbon-rich and 
silicate-rich dust. Lastly, we discovered a systematic change 
in the temperature profile for inner-most dust regions when 
the dust shell optical depth rises above T2.2 M m > 2. This 
observed lowering of the central dust temperatures could be 
naturally explained as a consequence of shadowing caused 
by clumpy dust formation on spatial scales smaller than our 
angular resolution, but other possibilities should be further 
explored as well. 
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Table 1. Properties of NIRC Camera Infrared Filters. Reference: The NIRC Manual. 



Name 


Center Wavelength 


Bandpass FWHM 


Fractional 




A (/im) 


AA (/an) 


Bandwidth 


Fell 


1.6471 


0.0176 


1.1% 


H 


1.6575 


0.333 


20% 


K 


2.2135 


0.427 


19% 


Kcont 


2.25965 


0.0531 


2.3% 


CH4 


2.269 


0.155 


6.8% 


PAHcs 


3.0825 


0.1007 


3.3% 



Table 2. Basic Properties of Targets 



Source 


RA (J2000) 


Dec (J2000) 


V 




H a 


K a 


Spectral 


Names 






mag 


mag 


mag 


mag 


Type 


AFGL 230 


01 33 51.21 


+62 26 53.5 




16.747 


11.232 


7.097 


mo 


AFGL 2019 


17 53 18.9 


-26 56 37 


20.2( 2 ) 


6.338 


4.035 


2.616 


M8W 


AFGL 2199 


18 35 46.48 


+05 35 46.5 




8.04 


4.85 


2.701 


M< 6 ) 


AFGL 2290 


18 58 30.02 


+06 42 57.7 




13.169 


8.966 


5.862 


M( 6 ) 


CIT 1 


00 06 52.94 


+43 05 00.0 


9.00W 


3.041 


1.829 


1.115 


M9W 


CIT 3 


01 06 25.98 


+12 35 53.0 




7.45 


4.641 


2.217 


M9W 


vl300 Aql 


20 10 27.87 


-06 16 13.6 


20« 


6.906 


3.923 


2.059 


M« 


AFGL 1922 


17 07 58.24 


-24 44 31.1 




12.244 


9.181 


6.342 


C( 3 ) 


AFGL 1977 


17 31 54.98 


+17 45 19.7 


9.9< 4 ) 


10.536 


7.994 


5.607 


CO) 


AFGL 2135 


18 22 34.50 


+27 06 30.2 




9.043 


6.002 


3.643 


C (i) 


AFGL 2232 


18 41 54.39 


+17 41 08.5 


9.7« 


5.742 


3.444 


1.744 


C (i) 


AFGL 2513 


20 09 14.22 


+31 25 44.0 




8.229 


5.705 


3.69 


c(i) 


AFGL 2686 


20 59 08.88 


+27 26 41.7 


20« 


9.112 


6.268 


4.075 


Ce(i) 


AFGL 4211 


15 11 41.89 


-48 20 01.3 




10.711 


7.751 


5.154 


0(3) 


IRAS 15148-4940 


15 18 22.05 


-49 51 04.6 


11.8« 


5.297 


3.071 


1.696 


c(i) 


IY Hya 


10 17 00.52 


-14 39 31.4 


14d) 


5.919 


3.666 


1.964 


C (5) 


LP And 


23 34 27.66 


+43 33 02.4 




9.623 


6.355 


3.859 


C (i) 


RV Aqr 


21 05 51.68 


-00 12 40.3 


11.5W 


4.046 


2.355 


1.239 


C (5) 


vl899 Cyg 


21 04 14.8 


+53 21 03 


15.6C 1 ) 


10.84 


8.693 


6.596 


C8( 5 ) 


V Cyg 


20 41 18.2702 


+48 08 28.835 


7 . 7 (i) 


3.096 


1.273 


0.117 


C(i) 



a These magnitudes (from 2MASS) are merely representative since the targets are variable. See Table 3 for our new photometry. 
Note: Horizontal line separa tes oxygen-rich (to p) from carbon-ric h (b ottom). 

References: (1) Simbad (2) |Monet[||l998r i. (3) iBuscombe] (|l998h . (4) lEgret etaD jl992j ). f5) ISkifj l|2007l ). (6) IVolld l|2002l ) (7) 
Garcia-Hcriian dez et al.l | |2007| ) 
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Table 3. Journal of observations and derived photometry 



Target 


Date (s ) 
(UT) 


Filter 


Aperture 
Mask 


Magnitude 


Calibrator 
Names 




AFGL 230 


1997 Dec 


k 

PAHcs 


FFA 

KL Relation* 


8.34±0. 1 
5. 11±0.2 


X Cas 






2002 Jul 


k 


FFA 


S.99±0.1 


HD 9878 








PAHcs 


FFA 


5.90±0.3 


HD 9329 




AFGL 2019 


2000 Jun 


CH4 


annulus 36 


2.48±0.1 


HD 163428 








h 


annulus 36 


3.84± 0.1 


HD 156992 








PAHcs 


annulus 36 


1.52±0.1 


HD 163428 




AFGL 2199 


1998 Apr 


CH4 


annulus 36 


2.99±0.1 


HD 170137 








PAHcs 


annulus 36 


1.80±0.1 


HD 170137 




AFGL 2290 


1998 Jun 


CH4 


annulus 36 


4.72±0.1 


HD 173074 








PAHcs 


annulus 36 


2.60±0.1 


HD 173074 






1999 Apr 




annulus 36 




HD 173833 








k 


annulus 36 


6.19±0.32 


HD 231437 








PAHcs 


annulus 36 


3.29±0. 1 


HD 173833 




CIT 1 


2000 Jun 


CH4 


annulus 36 


2 . 60 zb . 1 


A And 








h 


annulus 36 


4. 18±0 . 25 


HD 222499 








PAHcs 


annulus 36 


1.51±0. 1 


A And 




CIT 3 


1997 Dec 


Kcont 


annulus 36 


1 .08±0. 1 


5 Psc 








PAHcs 


annulus 36 


-0. 14±0. 1 


5 Psc 






1998 Sep 


CH4 


Golay 21 


2.45±0.1 


6 Psc 








PAHcs 


Golay 21 


1.04±0.1 


S Psc 




vl300 Aql 


1998 Jun 


CH4 


annulus 36 


1 .39±0. 1 


HD 189114 










annulus 36 


3.29±0.25 


HD 192464 








PAHcs 


annulus 36 


0.60±0. 1 


HD 189114 






1999 Jul 


kcont 


annulus 36 


2 . 02±0 . 1 


SAO 14382 








PAHcs 


annulus 36 


0.86±0. 1 


SAO 14382 




AFGL 1922 


2000 Jun 


k 

PAHcs 


annulus 36 
KL relation* 


6.34±0.25 
3.62±0.25 


HD 156992 






2001 Jun 


k 

PAHcs 


annulus 36 
KL relation* 


4.90±0.1 
2.32±0.25 


HD 158774 




AFGL 1977 


1998 Jun 


CH4 


annulus 36 


4.19±0.1 


HD 158227 








li 


annulus 36 


7.05±0.1 


HD 158227 








PAHcs 


annulus 36 


1.84±0.1 


HD 157049 






1999 Apr 


CH4 


annulus 36 


2.77±0.1 


HD 157049 








PAHcs 


annulus 36 


0.59±0.1 


HD 157049 




AFGL 2135 


2001 Jun 


k 


annulus 36 


3.29±0.1 


HD 168366, 


HD 181700 






PAHcs 


annulus 36 


1.27±0.3 


HD 177716 




AFGL 2232 


1998 Jun 


CH4 


annulus 36 


2.04±0. 1 


HD 158227 








h 


annulus 36 


4. 12±0. 1 


HD 158227 








PAHcs 


annulus 36 


0.68±0.1 


HD 157049 








CH4 


Golay 21 


2.28±0.3 


HD 168720 








PAHcs 


Golay 21 


0.94±0.3 


HD 168720 






1999 Apr 


CH4 


annulus 36 


1.06±0.1 


HD 173833 








PAHcs 


annulus 36 


-0.38±0.1 


HD 173833 




AFGL 2513 


1998 Sep 


h 


annulus 36 


6.58±0.1 


HD 196241 








CH4 


annulus 36 


4.03±0.1 


HD 200451 








PAHcs 


annulus 36 


3. 16±0.3 


e Cyg 






1999 Jul 


CH4 


annulus 36 


2.90±0. 1 


HD 188947 








PAHcs 


annulus 36 


1.70±0.1 


HD 188947 




AFGL 2686 


1998 Sep 


CH4 


annulus 36 


2.95±0.1 


HD 200451 








h 


annulus 36 


5. 82 ±0.21 


HD 200451 








PAHcs 


annulus 36 


1.01±0. 1 


e Cyg, A An 


d 




1999 Jul 


CH4 


annulus 36 


5. 13±0.1 


HD 188947 








PAHcs 


annulus 36 


2.92±0. 1 


HD 188947 








h 


annulus 36 


8.48±0.3 


HD 198330 




AFGL 4211 


2000 Jun 


CH4 


annulus 36 


3.62±0.3 


HD 137709 








PAHcs 


annulus 36 


1.42±0.3 


HD 137709 






2001 Jun 


k 

PAHcs 


annulus 36 
KL relation* 


4.70±0.1 
2.64±0.2 


HD 137709 




IRAS 15148-4940 


2001 Jun 


CH4 


annulus 36 


1.25±0.3 


HD 137709 








k 


annulus 36 


1.30±0.3 


HD 137709 








PAHcs 


annulus 36 


1 . 7 1 ± . 1 


HD 136422 




IY Hya 


1999 Apr 


CH4 


annulus 36 


2.08±0.1 


HD 87262 








PAHcs 


annulus 36 


1.37±0.1 


H Hya 




LP And 


1998 Sep 


CH4 


annulus 36 


3.89±0.1 


HD 222499, 


X And 






h 


annulus 36 


7.05±0.25 


HD 222499 








PAHcs 


annulus 36 


1.72±0.1 


A And 






1999 Jul 


CH4 


Golay 21 


4.01±0.1 


a Cas 








PAHcs 


Golay 21 


l.SOzbO. 1 


a Cas 






1999 Jan 


CH4 


Golay 21 


3.26±0. 1 


a Cas 








PAHcs 


Golay 21 


1. 18±0.1 






RV Aqr 


1999 Jul 


CH4 


Golay 21 


1.23±0.25 


SAO 143482 


3 Aqr 






PAHcs 


Golay 21 


0.56±0.25 


SAO 143482 


3 Aqr 




1998 Jun 


CH4 


Golay 21 


1.52±0.1 


HD 196321 








PAHcs 


Golay 21 


1.15±0.1 


HD 196321 




vl899 Cyg 


1998 Jun 


CH4 


annulus 36 


5.53±0.1 


HD 202897 








h 


annulus 36 


7.87±0.3 


HD 200817 








PAHcs 


annulus 36 


3. 71 ±0.1 


HD 202897 






1999 Jul 


k 

PAHcs 


annulus 36 
KL relation* 


G. id ±0.1 
4.72±0.2 


HD 198661 




V Cyg 


1998 Jun 


fcii 


annulus 36 


2.59±0.1 


HD 192909 








kcont 


annulus 36 


0.53±0.1 


HD 192909 








CH4 


Golay 21 


0.50±0.1 


HD 192909 








PAHcs 


annulus 36 


0.26±0.1 


HD 192909 








PAHcs 


Golay 21 


0.19±0.1 


HD 192909 






1999 Apr 


CH4 


Golay 21 


0.15±0.1 


5 Cyg 








PAHcs 


Golay 21 


-0.25±0.1 


£ Cyg 






2001 Jun 


CH4 


annulus 36 


-0.27±0.1 


£ Cyg 








PAHcs 


annulus 36 


-0.69±0.1 


£ Cyg 





* This point was extrapolated from another epoch for the same star and assigned an error of 0.2 mag. 
Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Table 4. Basic Properties of Calibrators 



Calibrator 


.1 




11 


K 




mag 


mag 


mag 


HD 168720 


1 


.79 


0.875 


0.870 


HD 170137 


3 


.476 


2.737 


2.230 


e Cyg 





.641 


0.2 


0.1 


HD 200451 


1 


101 


3.231 


2.840 


HD 231437 


5 


.027 


3.958 


3.693 


HD 173833 


3 


.488 


2.647 


2.1 


HD 158227 


5 


.626 


4.984 


4.812 


HD 157049 


1 


.975 


1.149 


.830 


HD 168366 


5 


.049 


4.535 


4.255 


HD 181700 


3 


.938 


2.993 


2.735 


SAO 143482 


1 


.665 


0.790 


0.573 


HD 189114 


3 


.212 


2.030 


1.953 


HD 137709 


2 


.232 


1.532 


1.331 


HD 222499 


i 


.641 


3.804 


3.627 


A And 


1 


.970 


1.4 


1.287 


HD 9878 


6 


.631 


6.730 


6.698 


HD 9329 


1 


.961 


4.381 


4.341 


HD 156992 


3 


.901 


3.123 


2.926 


HD 158774 


i 


.403 


3.451 


3.138 


HD 198611 


3 


.755 


2.862 


2.470 


HD 202987 


3 


.859 


3.067 


2.82 


3 Aqr 





.934 


-0.020 


-0.220 


HD 192909 


1 


.190 




0.180 


£ Cyg 


.995 


0.130 


-0.070 


HD 200817 


1 


.174 


3.721 


3.708 


HD 192464 


5 


.180 


4.176 


3.879 


a Cas 





.371 


-0.191 


-0.270 


H Hya 


1 


.216 


.506 


0.37 


HD 87262 


2 


.974 


2.052 


1.880 


HD 196321 


2 


.128 


1.361 


1.21 


HD 136422 








0.8 


S Psc 


2 


.031 


1.198 


.890 


HD 198330 


1 


.988 


4.159 


3.816 


HD 188947 


1 


.934 


1.438 


1.621 


X Cas 


3 


.019 


2.481 


2.311 


HD 163428 








1.6 


HD 196241 


4 


.19 


3.620 


3.090 



PAHcs 
mag 



Reference 



0.794 

2.16 

0.011 



2.02 
.684 

0.436 
1.908 
1.257 

1.245 

4.29 



2.75 
-0.338 
0.101 
-0.150 



-0.399 
0.28 

.98496 

0.535 

0.739 

1.561 

1.464 



'2003) 
'20031 
'20031 
'20031 
'20031 



Ncugebau<£_&_Jje^glrton (1969 
^Jejigebauer_ 
^Jeugebauer_ 



Leighton (,1969) 
Leighton (1969) 



Gullixson et al. (19831 
Gosnell et al. (19791 
extrapolation 



Price Murdock 



McWilliam fc Lambert (19841 . Cutri (20031 . Neugebauer fc Leighton 
Cutri (20031 . Neugebauer fc Leighton (19691 

Neugebauer fc Leighton (1969 1. Ghosh et al. (1984 1. Price fc Murdock 
Cutri 
Cutri 
Cutri 
Cutri 
Cutri _ 
Cutri (20031 
Cutri (20031 
Cutri (20031 
Cutri (20031 . 
Cutri (20031 .' 
Cutri (20031 
Johnson_e^il 
Cutri (20031 
Cutri (20031 . 
Cutri (20031 
Cutri (20031 
Cutri (20031 



19661 . Price fc Murdock (19991 . Selbv et al. (1988 



extrapolation 



J^awara^et^iL (1983) 
Cutri (2003), Ncugcbauer^j^Jje^ghton 
Cutri (20031 . Neugebauer fc Leighton (19691 
Carter (19901 

Johnson_e^il i (1966) , i Neugcbaucr &_Jjeighton (1969 ) , I^n£e_&_^irrdock 

Johnson et al. (19661 . Noguchi et al. (19811 
Cutri (20031 
Cutri (20031 

Voelcker (19751 . Alonso et al. (19941 



Price fe Murdock 



Price fc Murdock (19991 
Price fc Murdock (1999 
Price fc Murdock (1999 



Cutri (2003) 
Cutri (20031 
Cutri (20031 
Price (19681 
Cutri (20031 
Cutri (20031 

Noguchi et al. (19811 . Elias et al. (19821 . Glass 
Cutri (20031 . Neugebauer fc Leighton (1969 
White fc Wing (19781 . Humphreys fc Nev (19741 
Morel fc Magnenat (19781 . Cutri (2003) 



Jolmson_et i ^il i (1966) 
Nejigel2auer i ^Ijeighton (1969) 
Ncugcbaue^^Jje^ghton (1969) 
Eggcn (19691 
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Table 5. Results from Circularly-Symmetric Gaussian Models 

Target Datc(s) Filter 



^0 

(±0.05 



FWHM 

(mas) 



x'/DOF 



AFGL 230 


1997 


Dec 


k 


FFA 





71 


32±3 


0.23 




2002 


Jul 


k 


FFA 





54 


34±3 


0.34 








PAHcs 


F FA 





74 


33±2 


0.05 


AFGL 2019 


2000 


Jun 


CH4 


annulus 36 





96 


10 + ? n 
— 10 


0.31 








h 


annulus 36 





90 


9±4 


0.65 








PAHcs 


annulus 36 





Or, 


2 1 ± 3 


0.27 


AFGL 2199 


199S 


A pi' 




annulus 36 
















PAH 


annu us 






22 + 3 


n "i i- 


AKGL 2290* 


1998 


ub 




annu us , 






22 + 4 


'. ^ 








PAH 


annu us , 






27+3 


n ' r n 




1999 


A p c 




annul s 
























'J To 


' r 1 








PAHcs 


annulus 36 
annu us 


- 


83 


36±'> 


36 


C't 1* 


9000 


UB 




annu us 






1 r _|_ r 


, ' ^ 
















14 + 3 










PAHc 


annulus 36 




o'i 


■'•()~ J 


46 


Cit 3* 


1997 


D 


kcont ^ 


annulus 36 
annu us 




V) 




1 1 








PAH 


annu 8 


- 


. 




',. 




1998 


bep 




° ay 




'■■n 




n ' ~ i- 








PAH 


° & T 





0) 






vl300 Aql* 


1998 


UB 










7 , T ( " 












annulus ^ 




Vi 


1 4 ± 3 


11 








PAHcs 


annulus 36 


o 


8 1 


23 + 3 


. . r >o 




1999 


Jul 


kcont 


annulus 36 





87 


18 + 5 


0.39 








PAHcs 


annulus 36 





90 


21±3 


0.52 


AFGL 1922 


2000 


JUB 


k 


annulus 36 


o 


76 


24 + 4 


0.88 




2001 


JUB 


k 


annulus 36 





83 


29 + 4 


0.76 








PAHcs 


annulus 36 





95 


58 + 2 


0.43 


AFGL 1977* 


1998 


JUB 


CH4 


annulus 36 





78 


24 + 4 


0.26 








h 


annulus 36 





76 


17 + 3 


0.68 








PAHcs 


annulus 36 





94 


34 + 2 


0. 11 




1999 


Apr 


CH4 


annulus 36 





96 


29 + 4 


0.25 








PAHcs 


annulus 36 





89 


52 + 2 


0.26 


AFGL 2135 


2001 


Jun 


k 


annulus 36 





66 


17+5 


0.49 




2001 


JUB 


PAHcs 


annulus 36 





50 


34 + 2 


0.13 


AFGL 2232* 


1998 


Jun 


CH4 


annulus 36 





83 


18 + 5 


1.32 








h 


annulus 36 





81 


14+3 


0.50 








PAHcs 


annulus 36 





00 


33 + 2 


0.56 








CH4 


Golay 21 





91 


20 + 5 


0.17 








PAHcs 


Golay 21 





00 


34 + 2 


0.11 




1999 


Apr 


CH4 


annulus 36 





69 


44 + 3 


0.19 








PAHcs 


annulus 36 





86 


42 + 2 


0.12 


AFGL 2513* 


1998 


Sep 


h 


annulus 36 


1 


00 


— 1 


1.15 








CH4 


annulus 36 





94 


-. „ + 6 

10 -io 


0.18 








PAHcs 


annulus 36 


1 


00 


16±4 


0.70 




1999 


Jul 


CH4 


annulus 36 


1 


00 


11+6 
— 9 


0.32 










annulus 36 





'!r 


~ J * ' 


0.36 




998 


p 




























7 










PAHcs 


annulus 36 


.. 


■NO 


or 4_n 
35 + 2 


()..-, 6 




1999 


Jul 


A 4 


annulus 36 


n 






0.()8 








PAHcs 


annulus 36 






o«io 


0. 11 




























annulus 36 


;. 






' 








CS 


annu us 3 






o - i *- 






.i"" 1 


. Ull 


k 


annulus Jt) 











TR A O 1 ^ 1 ,1 Q A (1A 11 


00 






annulus 36 


1 




i ^i- 


nil 










annulus Jt) 







13 + 


0.50 








PAHcs 


annulus 36 




86 


25 + 3 


0.30 


IY Hya 


1999 


Apr 


a 4 


annulus 36 





88 


14 + 6 


0.32 










annulus 36 






Or T ^ 


0.28 






Sep 






















annulus 36 




tu . 




n r r 


















■ "i 1 




1999 


I ] 




Gala 21 


j.j 


89 


24 + 4 


I) 00 








PAHcs 


Golay 21 





79 


48 + 2 


0.50 




1999 


Jan 


CH4 


Golay 21 





70 


25 + 4 


2. 11 








PAHcs 


Golay 21 





66 


35 + 2 


0.45 


RV Aqr* 


1999 


Jul 


CH4 


Golay 21 


1 


00 


8+8 


o.i 6 








PAHcs 


Golay 21 





96 


26 + 3 


0.21 




1998 


Jun 


CH4 


Golay 21 





98 


12 + 8 


0.13 








PAHcs 


Golay 21 


1 


00 


27+3 


0.36 


vl899 Cyg 


1998 


Jun 


CH4 


annulus 36 





88 


18 + 5 


0.35 








h 


annulus 36 





86 


16 + 3 


0.39 








PAHcs 


annulus 36 





92 


22 + 3 


0.32 




1999 


Jul 


k 


annulus 36 





93 


15 + 5 


0.62 


V Cyg 


1998 


Jun 


feii 


annulus 36 





87 


14+3 


0.92 








kcont 


annulus 36 





96 


1 6 ± 5 


1.15 








PAHcs 


annulus 36 





90 


34+2 


0.51 








CH4 


Golay 21 


1 


00 


18 + 5 


0.35 








PAHcs 


Golay 21 





92 


38 + 2 


0.10 




1999 


Apr 


CH4 


Golay 21 





93 


17 + 5 


0.15 








PAHcs 


Golay 21 





86 


38 + 2 


0.12 




2001 


Jun 


CH4 


annulus 36 





82 


19 + 5 


0.26 








PAHcs 


annulus 36 





83 


42 + 2 


. 1 5 



* target is asymmetric, see Table 7 fo 
Note: Horizontal line separates oxygen 



further details 
rich (top) from 
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Table 6. Results from Central Point plus Circularly-Symmetric Gaussian Models 

Target Date(s) Filter Aperture f Point f Gauss FWHM x 2 /DOF 

(±.05) (±0.05) (mas) 



AFGL 230 


1997 


Dec 


k 


FFA 





24 





52 


17±7 


0.21 




2002 


.1 ul 


k 


FFA 





30 





12 


98±7 


0.15 








PAHcs 


FFA 





50 





32 


86±8 


0.42 


AFGL 2019 


2000 


.Tun 


CH4 




iulus 36 





86 





14 


51 + 3° 
— 45 


0.28 








h 


an 


aulus 36 





00 





83 


!±1 


0.83 








PAHcs 


an 


iulus 36 





45 





51 


31±5 


0.27 


AFGL 2199 


1998 


A P L 


CH4 


an 


iulus 36 





38 





5 1 


19±9 


0.22 








PAHcs 


an 


aulus 36 





36 





68 


30±4 


0.43 


AFGL 2290* 


1998 


Jun 


CH4 




1U | US ^ 





44 





38 


46±12 


0.26 








PAHcs 


an 







55 





38 


68±7 


0.61 




1 9 9 9 


A p i' 


CH4 




m us 


Q 


31 


o 


51 


6 6 — 9 


0.21 








k 









34 


o 


56 


66±9 


0.31 








PAHcs 




1U l 





25 





60 


49±3 


0.35 


CIT 1* 


2000 


.Tun 


CH4 




1U ; U 8 





00 


o 


92 


14±6 


0.35 








h 




m us 





47 


o 


19 


23±6 


0.36 








PAHcs 









00 


o 


94 


20±4 


0.46 


CIT 3 * 


1997 


Dec 


kcont 


an 







58 


o 


50 


53 = 13 


0.23 








PAHcs 




iulus 36 





36 





62 


60±4 


0.02 




1998 


Sep 


CH4 


c n 


y & 2 1 





50 


o 


44 


40±10 


0.03 








PAHcs 


Go 


lay 21 





46 





47 


50±5 


0.20 


vl300 Aql* 


1998 


.Inn 


CH4 




aulus 36 





64 





23 


43±20 


0.40 








h 


an 


aulus 36 





45 





38 


25±7 


0.39 








PAHcs 


an 


aulus 36 





64 





34 


80±10 


0.30 




1999 


.Jul 


kcont 


an 


aulus 36 





00 





87 


18±5 


0.39 








PAHcs 




aulus 36 





17 





73 


24±4 


0.52 


AFGL 1922 


2000 


Jun 


k 


an 


3U . US 





12 





39 


47 = 11 


0.83 
































PAHcs 






o 


51 


o 


17 


105±6 


0.41 


AFGL 1977* 


1998 


.Inn 


CH4 


an 


au us 36 





36 





45 


1 1 = 9 


0.22 








h 


an 


aulus 36 





5 





33 


43 -13 


0.59 








PAHcs 


an 


aulus 36 





42 





58 


56±4 


0.36 




1999 


Apr 


CH4 









34 





69 


43±7 


0.17 








PAHcs 


■ui 


aulus 36 





2 5 





74 


45±3 


0.17 


AFGL 2135 


2001 


Jun 




an 


aulus 36 


Q 


00 





66 


17=5 


0.49 




2001 


.Tun 


PAHcs 


an 


aulus 36 


o 


28 


o 


27 


78±6 


0.10 


AFGL 2232* 


1998 


Jun 


CH4 


an 


aulus 36 





56 





33 


1-1=1-1 


0.47 








h 


an 


aulus 36 





00 





81 


14±3 


1.32 








PAHcs 


an 


aulus 36 





50 





49 


66±6 


0.49 








CH4 


Go 


lay 21 





52 





45 


39±10 


0.10 








PAHcs 


Go 


lay 21 





38 





60 


51±4 


0.09 




1999 


Apr 


CH4 


an 


aulus 36 





23 





56 


72±7 


0.08 








PAHcs 


an 


aulus 36 





26 





64 


58±3 


0.07 


AFGL 2513* 


1998 


Sep 


h 


an 


aulus 36 





03 


1 


00 


1 + 9 


1.09 








CH4 




m : us 


o 


81 


o 


14 


39±39 


0.17 








PAHcs 


an 


au us 


o 


03 


1 


00 


10—1 


0.68 




1999 


Jul 


CH4 


an 


aulus 36 





00 





92 


,+13 
1 -1 


0.65 








PAHcs 


an 


aulus 36 





67 





37 


65±9 


0.24 


AFGL 2686 


1998 


Sep 


CH4 




1U [ US 





30 





63 


42±7 


0.39 








h 


an 


1U . US 





24 





69 


35±4 


0.43 








PAHcs 






Q 


39 





5 6 


58±4 


0.29 




1999 


.lul 


CH4 


an 


1U . US 


o 


21 


o 


73 


32±5 


0.67 








PAHcs 









29 





64 


45±3 


0.42 










an 


aulus 36 





6 3 





3 


,,,+ 15 


0.77 


AFGL 4211 


2000 


Jun 


CH4 


an 


aulus 36 





36 





54 


59±10 


0.27 








PAHcs 


an 


aulus 36 





16 





70 


S9±4 


0.06 




2001 


."hm 


k 


an 


aulus 36 





47 





28 


86 ±19 


0.45 


IRAS 15148-4940 


2001 


.1 1111 


CH4 


an 


aulus 36 


o 


62 


Q 


18 


44=23 


0.39 










an 


aulus 36 


o 


52 


o 


31 


25 = 13 


0.59 








PAHcs 


an 


aulus 36 


Q 


52 





38 


49 ±7 


0.37 


IY Hya 


1999 


A p i 


CH4 


an 


aulus 36 


Q 


00 


o 


88 


14=6 


0.32 












aulus 36 














LP And* 


1 9 9 S 


Se 

P 


CH4 


an 


aulus 36 


o 


■12 


o 


19 


49= 10 


0.41 












aulus 36 


o 


37 


o 


34 


39 = 9 


0.52 








PAHcs 




aulus 36 


Q 


3 





65 


74=4 


0.35 




1999 


Jul 


CH4 


Go 


lay 21 





47 





49 


47±11 


0.91 








PAHcs 


Go 


lay 21 





38 





51 


85±5 


0.42 




1999 


Jan 


CH4 


Go 


lay 21 





47 





49 


47±11 


0.91 








PAHcs 


Go 


lay 21 





35 





51 


85±5 


0.42 


R.V Aqr* 


1999 


.Jul 


CH4 


Go 


lay 21 





00 





96 


n +13 
i -l 


0.30 








PAHcs 


Go 


lay 21 





33 





64 


34±4 


0.21 




1998 


.7un 


CH4 


Go 


lay 21 





74 





25 


30±20 


0.12 








PAHcs 


Go 


lay 21 





43 





62 


42±4 


0.33 


vl899 Cyg 


1998 


Jun 


CH4 




aulus 36 





32 





57 


24±7 


0.35 








h 




aulus 36 





72 





23 


75±15 


0.30 








PAHcs 




aulus 36 





65 





23 


57±13 


0.31 




1999 


Jul 


k 




aulus 36 





76 





23 


52±24 


0.59 


V Cyg 


1998 


Jun 


feii 




aulus 36 





31 





56 


1S±5 


0.91 








kcont 




aulus 36 





51 





17 


26±9 


1.14 








PAHcs 




aulus 36 





00 





90 


34±2 


0.51 








CH4 


G. 


lay 21 





58 





17 


35±10 


0.30 








PAHcs 


Go 


lay 21 





28 





67 


52±3 


0.06 




1999 


A P L 


CH4 


Go 


lay 21 





52 





43 


30±9 


0.14 








PAHcs 


Go 


lay 21 





31 





60 


57±4 


0.07 




2001 


.Inn 


CH4 




aulus 36 





50 





36 


40±11 


. 2 1 








PAHcs 




aulus 36 





28 





61 


63±4 


0.09 



* target is asymmetric, sec Tabic 7 for further details 

Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Table 7. Results from 2-dimcnsional Gaussian Models 



Target 


Date(s) 


Filter 


FWHM major 
FWHM minor 


<PA> 


AFGL 2290 


1998 Jun 


CH4 
PAHcs 


1.23 
1.24 


58±20 




1999 Apr 


CH4 
k 

PAHcs 


1.24 
1.24 
1.06 




CIT 1 


2000 Jun 


CH4 
h 

PAHcs 


1.13 
1.14 
1.11 


133±3 


CIT 3 


1997 Dec 


Kcont 
PAHcs 


1.19 
1.06 


151±9 




1998 Sep 


CH4 
PAHcs 


1.28 
1.04 




V1300 Aql 


1998 Jun 


CH4 
h 

PAHcs 


1.34 
1.14 
1.19 


108±13 




1999 Jul 


kcont 
PAHcs 


1.34 
1.31 




AFGL 1977 


1998 Jun 


CH4 
PAHcs 


1.11 
1.31 
1.06 


71±18 




1 qqq Arir 


CH4 
PAHcs 


1.10 
1.21 




AFGL 2232 


1998 Jun 


CH4 

PAHcs 

CH4 

PAHcs 


1.37 

l.oo 

1.19 
1.10 
1.08 


94±10 




1999 Apr 


CH4 
PAHcs 


1.22 
1.05 




AFGL 2513 


1998 Sep 


h 


unresolved 








CH4 


1 5 


61±8 






PAHcs 


1.39 






1 qqq Tnl 


CH4 
PAHcs 


1.38 
1.2 




LP And 


1998 Sep 


CH4 
h 

PAHcs 


1.46 
2.03 
1.39 


108±6 




1999 Jul 


CH4 
PAHcs 


1.64 
1.36 






1999 Jan 


CH4 
PAHcs 


1.82 
1.20 




RV Aqr 


1999 Jul 


CH4 
PAHcs 


1.49 
1.23 


122±24 




1998 Jun 


CH4 
PAHcs 


1.35 
1.22 





* Sources missing from this list were found to have circularly-symmetric dust shells (within errors). Note: PA is the mean position 
angle of the major axis (degrees East of North) for all filters and epochs. 
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Table 8. Results from DUSTY Radiative Transfer Model 



Target 


Date(s) 


Tdust 


T2.2/Z m 


R* 


^(1 fcpc) 


X 2 /DOF 








(K) 




(mas) 


(103 Lq) 




AFGL 230 


1997 


Dec 


800Ig 


a n+0.9 

4 - y 18.7 


1 e: + 0.5 

1-516.3 




0.26 




2002 


Jul 


^+400 
04U_ 110 


7 4+1-6 
'• 4 -1.2 


4 1+ 4 - 5 
4 - 1 -2.1 


o | +108 
ol_24 


O.OO 


AFGL 2019 


2000 


Jun 


ilyu -250 


Q9+0.23 
u - yz -0.12 


q pr + 0.7 
• 3 -°-0.3 


24lf 


0.54 


AFGL 2199 


1998 


Apr 


1130t 3 3 I° 


1 6+1-2 

1 - D -0.7 


q q+2.0 


21+f 


0.06 


AFGL 2290 


1998 


Jun 


B50tJS° 


q K+0.5 
°-°-0.5 


o 7 +0.8 
' -0.7 


26+1 2 


0.33 




1999 


Apr 




4.618:5 

i O+0.5 
l- z -0.2 


q Q+1.7 


29+31 


2.63 


CIT 1 


2000 


Jun 


iiau_ 230 


q c+0.8 


241J 2 


0.60 


CIT 3 


1997 


Dec 


1110+ 230 

iliu -140 


i 4+ - 7 


7 -»-0.6 


"6tiS 


0.34 




1998 


Sep 


1020+?™ 


1 q+0-7 


= n+0.8 
°- u -0.6 


48111 


0.29 


V1300 Aql 


1998 


Jun 


1080 +3 «! 


92+ ' 46 
u - yz -0.12 


K. o + l.O 
b -°-0.6 


64l 2 2 


0.50 




1999 


Jul 


1160l 34 8* 


i-60lg; 7 * 


tr 1+2.0 
°- L -0.8 


49l« 


0.11 



AFGL 1922 


2000 


Jun 


oc . n + 200 

850-60 


5 


o+0.7 

■ 6 -0.7 


1 


E+l-l 

°-0.8 


0-7+21 


1.44 




2001 


Jun 


850+1™ 


3 


q+0.5 

5 


5 


9+1.2 
Z -1.5 


51±1S 


0.39 


AFGL 1977 


1998 


Jun 


910l 8 9 ° 


2 


o + 0.'l 

°-0.2 


4. 


0+0.2 
U -0.6 


31+ 4 


1.69 




1999 


Apr 


990l 9 6 ° 


2 


c+0.5 
°-0.2 


6 


0+0.8 

A +65.7 
'+-4.6 


68I1 9 


0.65 


AFGL 2135 


2001 


Jun 


740! 3 ™ 


3 


9+4.2 

' z -1.4 


9 


167+1°,™° 


5.04 


AFGL 2232 


1998 


Jun 


lliollfo 


1 


6+ - 2 
D -0.2 


5 


1+0.5 
1-0.5 


48llo 


0.32 




1999 


Apr 


13001?™ 


2 


0+1.4 

'°-1.2 


9 


4 +5.0 


165+2f 


3.29 


AFGL 2513 


1998 


Sep 


1500l° 50 


1 


O+0.2 
y -0.2 


1 


7+0.9 
'-0.2 


= Q+7.3 
o -°-0.9 


0.57 




1999 


Jul 


1110+ 400 
H iu -200 


1 


9+0.9 
z -0.5 


3 


i+l.O 
1-0.4 


18ll 4 


0.30 


AFGL 2686 


1998 


Sep 


H10+ 140 
iliu -140 


2 


O+0.5 
'8-0.5 


5 


o + l.l 
°-0.8 


53l? 4 


1.85 




1999 


Jul 


8201™ 


3 


9+0.2 
z -0.2 


3 


-1+0.4 
1-0.3 


1911 


1.02 


AFGL 4211 


2000 


Jun 


880+™ 


3 


7+0.7 
'-0.S 


6 


7+0.3 
'-0.3 


85l 7 


0.73 




2001 


Jun 


8501^° 


4 


9+0.9 
z -0.7 


6 


-1+3.7 
1-2.0 


71+113 

'1-39 


3.20 


IRAS 15148-4940 


2001 


Jun 


9A0 + \f * 





9O+0.46 


4. 


fi +0.2 
°-0.9 


40lt 4 


2.47 


IY Hya 


1999 


Apr 


9601$ 





4fi +0.46 
' 40 -0.12 


4. 


9+0.1 
z -0.2 


33+2 


0.25 


LP And 


1998 


Sep 


8801™ 


3 


0+0.5 
u -0.2 


4. 


n+0.9 
' 8 -0.5 


13+j u 


1.49 




1999 


Jul 


8201™ 


3 


0+0.2 


4. 


0+O.6 
°-0.7 


44lH 


1.00 




1999 


Jan 


8801™ 


3 


9+0.5 
z -0.2 


6 


7+1.0 
' -0.7 


85i? 7 


1.67 


RV Aqr 


1999 


Jul 


1500l° 80 





46+0-23 
' 40 -0.23 


5 


4+0.7 
' 4 -0.4 


55+1 5 


1.09 




1998 


Jun 


1190l?i° 





9O+0.46 

••"-0.12 


5 


1+0.1 
1-0.7 


48l 2 3 


0.23 


vl899 Cyg 
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* This star has two regions which meet our 1-cr criteria for a best fit. The particular values shown were chosen for consistency, see the 
appropriate figure for more details. 

Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Figure 1. Best fit plots for AFGL230. The first row are figures for 
the epoch Dec97 and the second row is for Jul02. The first panel 
in each row shows a fit to the SED with our new photometry in- 
cluded with errors (2MASS points are plotted as squares in each 
frame for reference). The dashed line represents the contribution 
from the star, the dotted line represents dust contribution, the 
dash-dotted line represents the contribution from scattered light, 
and the solid line is the total flux. The second panel shows our 
DUSTY fits to the visibility data for each wavelength of obser- 
vations. The third panel shows the \ 2 /DOF surface, with bright 
areas showing the best-fitting region. The black contour denotes 
the 1-{<t} error . 



Figure 4. Best fit plots for AFGL2290. See Fig.l caption. 
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Figure 5. Best fit plots for CIT 1. See Fig.l caption. 




Figure 2. Best fit plots for AFGL2019. See Fig.l caption. 





Figure 6. Best fit plots for CIT 3. See Fig.l caption. 



Figure 3. Best fit plots for AFGL2199. See Fig.l caption. 
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Figure 7. Best fit plots for vl300 Aql. See Fig.l caption. For 
the 1999 Jul. epoch we chose the lower-right region as the best 
fit region because it is consistent with the best fit region for the 
1998 Jun. epoch. 
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Figure 8. Best fit plots for AFGL1922. See Fig.l caption. 




Figure 9. Best fit plots for AFGL1977. See Fig.l caption. 
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Figure 10. Best fit plots for AFGL2135. See Fig.l caption. 
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Figure 11. Best fit plots for AFGL2232. See Fig.l caption. 
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Figure 12. Best fit plots for AFGL2513. See Fig.l caption. 

10 13 [A^GL26&6 ' 




Figure 13. Best fit plots for AFGL2686. See Fig.l caption. 
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Figure 14. Best fit plots for AFGL4211. See Fig.l caption. 






Figure 17. Best fit plots for LP And. See Fig.l caption. 



Figure 15. Best fit plots for IRAS15148-4940. See Fig.l caption. 
We chose the lower-right region as the best fit region because in 
all other cases of multiple good fitting regions the one at low tau 
and high dust temperature was the consistent region. 






Figure 16. Best fit plots for IY Hya. See Fig.l caption. 




Figure 18. Best fit plots for RV Aqr. See Fig.l caption. 





Figure 19. Best fit plots for vl899 Cyg. See Fig.l caption. 
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Figure 20. Best fit plots for V Cyg. See Fig.l caption. For 
the 1998 Jun. epoch the best fitting region was chosen to be the 
lower-right region because it is consistent with the other epochs. 
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Figure 21. A plot of best fit T2.2^m versus the temperature 
at the inner edge of dust shell T<j ust . Open symbols are used for 
oxygen-rich dust shells and closed symbols are used for carbon- 
rich dust shells. 
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